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The prediction of performance of photovoltaic technologies is crucial, not 
only to improve the reliability and durability of these technologies but also to 
increase the confidence of investors and consumers in them. The accurate 
calculation of the degradation rate Dr (%) in real operating conditions under 
specific climatic stresses is, therefore, paramount. The present study provides 
a comparison of performance losses of two silicon PV technologies installed 
on the rooftop of the Higher School of Technology in La&youne-Morocco. 
The two systems are a polycristalline array (pc-Si: 1.82 kWp) and an 
amorphous array (a-Si: 1.55 KWp), which are grid connected. In the light of 
related performance gathered over three-year, the degradation rates of the two 
systems were estimated using four statistical methods under the open-source 
software R. The techniques engaged in this paper are: classical seasonal 
decomposition (CSD), holt-winters (HW), autoregressive integrated moving 
average (ARIMA), and seasonal and trend decomposition by LOESS (STL). 
The results obtained using those methods show that DR(%) varies between 
0.39% and 0.99% for pc-Si and between 0.29% and 0.64% for a-Si. The 
analysis of degradation accuracy shows that STL and CSD techniques provide 
results with high accuracy than ARIMA and HW for the two systems. The 
present study adds to knowledge on PV degradation under the subtropical 
desert climate of Ladyoune. 
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NOMENCLATURE 


a-Si: Amorphous Silicon 


ARIMA: Autoregressive integrated moving average 


c-Si: Crystalline Silicon 
CdTe: Cadmium telluride 


P*u: Effective peak power (W) 

Po: PV rated power (Wp) 

P* exp: Experimental peak power (W) 
pe-Si: Polycrystalline Silicon 


CIS: Copper indium-selenide 

CSD: Classical seasonal decomposition 

ES: Exponential smoothing 

Go: Solar irradiance under standard test conditions (W/m) 
He: Total solar radiation arriving at the surface of PV panels 
(Wh/m?) 

HIT: Hetero-junction with intrinsic thin-layer 

HW: Holt-winters 


Journal homepage: http://ijeecs.iaescore.com 


PID: Potential-induced degradation 

PR: Performance ratio (%) 

PRoc: DC- performance ratio (%) 

PRt: Temperature corrected performance ratio (%) 
Pt: Temperature-corrected power (W) 

PV: Photovoltaic 

PVUSA: Photovoltaics for utility scale applications 
Ro: Degradation rates (%) 


1234 O ISSN: 2502-4752 


I: Current (A) STC: Standard test conditions 

LID: Light-induced degradation Ta: Ambient temperature (°C) 

LLS: Linear least squares Tcell: temperature of the PV cell (°C) 
LOESS: Locally weighted scatterplot smoothing V: Voltage (V) 

LR: Linear regression Ya: Array yield (kWh/kWp) 

mce-Si: Monocristalline silicon Yf: Final yield (kWh/kWp) 

MJ-Si: Multi-junction silicon YOY: Year-on-year 

P: Power output (W) STL: Seasonal-trend decomposition by Loess 


KToe: Kilotonne of oil equivalent 


1. INTRODUCTION 

The need of energy in the world rises every year. According to the International Energy Agency (IEA), 
the total energy consumption in the world in 2018 was 9 937 703 KToe. Giving the increasing prices of fossil 
energies and dramatic environmental issues brought about them, renewable energies stand as one of promising 
alternatives to ensure energy supplies along with environment protection. As any other technology; to entice 
more investors and consumers; the cost, the reliability and the lifetime must be attractive. Therefore, many 
studies were established in order to lower the costs and to guarantee high reliability, efficiency and durability 
of photovoltaic (PV) systems [1]-[4]. 

The loss of performance and degradation of PV systems are generated by external as well as internal 
factors. External ones are generally: temperature cycles [5], humidity, solar irradiation, mechanical stress, 
soiling, and snow. As for internal ones, the main cause is the imperfect production. Those factors result in 
different type of failure modes who are responsible not only for power loss but also for safety issue [6]. 
The prediction of the degradation rate Dp (%) of PV systems with high accuracy, during outdoor operation and 
under different climate conditions is still incompletely explored. The intention of this work is to enrich the pool 
of knowledge on PV degradation: using a case study of two silicon PV technology systems operating outside 
under climate conditions of Laaéyoune City, in order to assess the performance degradation and determine which 
system is more suitable in such conditions. 

Performance loss springs up at different level: cell, module, array and system with various factors and 
degradation mechanisms at each level [7]. The paper exhibited different failure modes and showed their related 
impacts on the performance of PV module, as well as defining the risk priority number “RPN” for each defect 
who describes which one affect the system performance more [8]. This study showed also that hot-spot and 
encapsulant delamination, back sheet delamination, corrosion in solder bond fatigue are crucial failures 
connected with the security issue with lower RPN number <S50 [8]. In the literature, many methods are used in 
order to determine the performance loss of PV systems, which can be classified into two types: analytical 
models and statistical methods. The first ones are generated from the physical/chemical theories of a specific 
failure, by representing the mechanism involved in complex physical/chemical processes [9]. Statistical 
methods generally used are namely: simple linear regression (SLR) classical seasonal decomposition (CSD) 
holt-winters (HW) auto-regressive integrated moving average (ARIMA) and locally weighted scatterplot 
smoothing method (LOESS). The various choices of the measuring equipment, the data qualification, filtering 
criteria, the performance metric and the statistical method of estimation of the trend, spin off different results 
with varying uncertainty [7]. 

The estimation of the degradation rate Dr (%) of peak power of two grid-connected PV systems of 
two different technology having the same capacity of 24 kWp was conducted by using linear least square fitting 
method [10]. Dag and Buker [10] reveled that the thin-film technology (HiT) shows lower degradation rate 
with about 0.1% than the poly-crystalline (pc-Si) based technology within the range of 0.67% to 0.83%, 
respectively, after two and half year out-dour exposure under the semi-arid weather of Anatolia-Turkey. 
Another study from Solfs-Aleman et al. [11] was interested in exploring the linearity of the degradation rate of 
various thin film technologies in Jaén, Southern Spain. In this paper two methods were applied: classical 
seasonal decomposition (CSD) and year-on year (YOY) in order to discover the path of degradation of two PV 
field PVGCS #1(a-Si) and PVGCS #2 (a-Si/uc-Si) of 60 Wp and 110 Wp capacity, respectively. Four PV 
modules of four technologies were also studied -a-Si, a-Si:H/ wtc-Si:H, cadmium telluride (CdTe) and copper 
indium gallium selenide (CIGS)- under outside conditions. According to the first method, the degradation rates 
were around 1.32%, 1.22%, 1.50% for a-Si, a-Si/uc-Si fields and a-Si module respectively [11]. The a-Si/uc- 
Si PV module shows a lower power loss (~0.6%/year). The CdTe PV module, likewise, has a degradation rate 
who is nearly negligible according to CSD, while the YOY indicates that such parameter lies nearby 0.7%/year. 
Whereas the CIGS PV module exhibits an interesting improvement at a rate around—0.92%/year, as presented 
by the same study [11]. An example for utility scale PV system was explored by the study [12], in which, a 
system of 1MWp installed in a tropical semi-arid climate in Telangana State- India was investigated by 
applying LLS, CSD, HW and STL to the time series of performance ratio collected over four years of operation, 
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the degradation rate found was respectively 0.27 + 0.22%/year, 0.32 + 0.11%/year, 0.50 + 0.22%/year, and 
0.27 + 0.15%/year, according to the study. 

Other studies were interested in analyzing the degradation rate caused by specific failures. The paper 
reveals that after three years of outdoor exposure, under semi-arid climate in Bengurir-Morocco; the power 
degradation rate reached 4.57%/year in the modules with breakage and cracks [13]. After evaluating 144 
failure-survey-data collected from 18 countries, the most dominant defects encountered according to study [14] 
were: cell crack, potential induced degradation by shunts, defect bypass diodes and discoloring of the pottant. 
The mean degradation rates in moderate climate for cell crack, PIDs and discoloring of pottant were 5%, 16%, 
and 1% respectively. Another study underlines the main PV defects occurring after a longer-term operation of 
22 years under the prevailing weather in Seville-Spain. Lillo-Sanchez et al. [15] found that the important 
aspects of degradation occurred were: severe browning, milky pattern and oxidation of the metallization grid. 
That results in a degradation rate of the mean peak power of 1.4%/years. An attempt to enhance, more globally, 
the understanding of PV performance and degradation was made by the study [16], by providing a map of 
degradation mechanisms and degradation rate for mono-crystalline PV module. The map can be considered a 
guide for localizing zones with high climate stress. However, the development of a global map of the 
degradation for each PV technology under diverse operating and climatic conditions needs further 
contributions. This paper provides a case study of two grid-connected PV technologies operating outdoors 
under subtropical climate conditions. 

In light of this information, the calculation of the degradation rate DR (%) is important to lower the 
performance uncertainty and predict their evolution with high accuracy. Firstly, this paper will present the two 
systems under consideration as well as the four statistical techniques used (section 2). Secondly, this study will 
evaluate the performance loss rate of the two systems using the statistical methods (Section 3) and provide a 
comparison between the results given by these techniques and their accuracy. Finally, A comparison of results 
in this work with previous papers will also be mentioned. The present work is dedicated to exploring how the 
two PV systems perform under the prevailing weather in Ladyoune City, and how their performances degrade 
using different statistical techniques under R software. 


2. MATERIALS AND METHOD 
2.1. Materials description 

The two PV systems under the scope of the study, were installed in March 2018 as part of the 
“Propre.ma” project; who aims to develop a map of the productivity of different PV technologies over the 
Moroccan Kingdom; in the north-west of the ESTL, University IBN Zohr (latitude 27°07'50.4 "north and 
longitude 13°08'18.1" west) under semi-arid climate conditions. The installation is fundamentally formed by 
three PV technologies Poly-crystalline p-Si (1.82 kWp) and Mono-crystalline m-Si (2 kWp) and Amorphous 
silicon thin film technology a-Si (1.55 kWp) who are linked together in order to form a grid connected system 
of total power of 5.365kWp. A previous work was interested in determining the performance of these two 
silicon technologies in addition to mc-Si array [17]. In this study, only the pc-Si and a-Si will be considered. 

The amorphous PV system comprises ten amorphous silicon thin-film panels linked in series of 
155Wp each. As for the poly-crystalline PV system, it contains seven panels, each with power of 260 W linked 
also in series. The Table 1 below represents the technical specification of the two types of modules, as shown 
in the panel’s nameplate of each technology: 


Table 1. Technical specification of PV modules 


Trade mark Nexpower Solarworld 
Model NT-155AF Sunmodule plus SW 260 poly 
Solar cell Thin film, Amorphous Si Polycrystalline 
Maximum power At STC (Pm) 155 W 260 W 
Maximum power point voltage (Vmp) 65.9 V 31.40 V 
Maximum power current (Imp) 2.43 A 837A 
Open circuit voltage (Voc) 85.5 V 34.40 V 
Short circuit current (Isc) 2.56A 8.94 A 
Length 1.4m 1.675 m 
Width lm lm 
Weight 19.5 Kg 18Kg 


2.2. Performance metrics 

Performance metrics are generally used in order to readily compare different PV systems in different 
climates conditions. In this paper, we have based the estimation of performance loss Dr (%) on normalized 
performance metric, which is array yield (Y,). Monthly data are recorded over a monitoring period of three 
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years from January 2018 to December 2020, in order to establish time series of Ya. The expressions of 
performance parameters are based on definitions standardized by IEC 61724 [18]. 


2.2.1. Reference Yield (Yr) 

The reference yield Yr reflects the number of hours during which the solar radiation would need to be 
at reference irradiance levels to contribute the same incident energy. The array yield Yr has unit of hour (h). 
This yield can be calculated by dividing the total solar in-plane irradiation Ht (KWh/m2) by the reference in- 
plane irradiance GO, (1 KW/m2), as given by the following (1). 


Ht 

aor (1) 
2.2.2. Array yield (Ya) 

The array yield Ya is the daily array energy output Ea.«per kW of installed PV array Po. This yield 
represents the number of hours per day that the array would need to operate at its rated output power Po to 
contribute the same daily array energy to the system as was monitored [12]. It is calculated by using the 
expression: 


Y=" (2) 


2.3. Methodology of Statistical methods 

In order to evaluate the annual degradation rate of the two silicon PV technology under study (pc-Si 
and a-Si), four statistical methods are applied to the monthly time series of Array Yield Y, of each technology. 
The classical seasonal decomposition (CSD), holt winters exponential smoothing (HW), autoregressive 
integrated moving average (ARIMA) and seasonal and trend decomposition using LOESS (STL) techniques 
were implemented to calculate the degradation rate (Dp in units of %/ year) from the trend of PV module 
performance time series, using the slope of such a trend given by applying linear regression [7]. The simulation 
was done using R, an open-source Software. The annual degradation rate of PV grid-connected system can be 
given by (3): where a is the slope of the trend. 


Da(%) =a*12 (3) 


2.3.1. Classical seasonal decomposition (CSD) 

In the CSD method, the time series is decomposed to three elements: the trend, seasonal and remainder 
random element. The trend is a result of the application of a moving average centered on two steps. Giving a 
moving average of 2k, where k is the order of the moving average, the centered average at time t which 
represents the trend (Tt) at time t, which is given by (4) [7]: 


1/1 t+m-1 1 t+m 
pa apes X) 
2 k i=t—m k i=t—m+1 (4) 
where: m=k/2 is defined as half the width of a moving average [9]. 
As for the seasonal element, it is given by subtracting the trend from the initial time series data, 


according to the (5) [9]: 
Se=¥%; -T, (5) 


the degradation rate Dr (%) is, then, calculated from the trend element by using linear regression. The Dy (%) 
is calculated by multiplying the slope ‘a’ of the linear regression curve of the extracted trend by 12 as presented 


in (3) [9]. 


2.3.2. Holt-winters model (HW) 

The Holt-winters (HW) technique is derived from a triple exponential smoothing method. This model- 
based method is applied to the time series and takes into account seasonal variations in addition to trends. The 
model establishes the trend, seasonal and residuals components by the application of centered moving average, 
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through the minimization of the squared one-step ahead error [7]. The HW method for PV systems follows the 
general additive model as shown in (6) [9], [19]: 


Ye+at = Ly + by + Stes (6) 


where I, is the level component, b, the slope component and S,_,,1 the relevant seasonal component, s the 
seasonal period, their expressions are given in (7)-(9) [9]: 


ly = AQ: — Sts) + 1 — A)Ce-a + D¢-1) (7) 
by = BCL; — ly-1) + (1 — B)by-4 (8) 
Se = CMe — Le-a — Dp) + A - C)Sy-5 (9) 


where A, B and C values are within [0,1] interval. 


2.3.3. Autoregressive integrated moving average (ARIMA) 

ARIMA is a model, which holds various methods in a multiplicative way, and can be abbreviated as 
ARIMA (p, d, q) (P, D, Q). Where p is the auto-regressive, d the differencing, and q the moving average order. 
As for P, D, and Q, they represent the seasonal autoregressive, the seasonal differencing, and the seasonal 
moving average order, respectively. [9] One of the most important advantages of the model is its flexibility to 
deal with seasonal fluctuations, errors, outliers, and level shifts, in a suitable way [9]. The general model of 
ARIMA can be given by the following expression (10). 


PT) bs(T*)V VE y¥¢ = O(T)O5(T Jee (10) 


Here T represents the delay operator, (T) = (1 = $1T—..—pT” )represents an autoregressive 
polynomial in T of degree p, @(T'*) is an autoregressive polynomial in Ts of degree PS. @(T) is a moving 
average polynomial in T of degree q, 0,(T*) is a moving average polynomial of degree Qs in Ts. The operator 
V4= (1 —T)* is a non-seasonal differencing operator, and V? is a seasonal differencing operator and captures 
non-stationarity in the corresponding locations in consecutive periods [9]. 


2.3.4. Seasonal-trend decomposition using LOESS (STL) 

STL decomposition is a non-parametric method, which estimates non-linear relationship by using a 
locally weighted regression instead of centered moving average to extract the trend. Thus, the method is more 
robust and less influenced by missing values and outliers [9]. Furthermore, it is flexible in defining the amounts 
of fluctuations in the trend and seasonal components. In addition, it has the flexibility in specifying the number 
of observations per cycle of seasonal component to any integer greater than one [20]. As CDS method, STL 
model aims to decompose the time series to three components as define in (11). 


Y, = T; +5, 4+ Rt t=1,2 ecudnareceies ~n (11) 


Where T;, S; and R; are the trend, the seasonal and the residual components, respectively [9]. 


3. RESULTS AND DISCUSSION 
3.1. Variation and seasonality of parameter metric 

Figure | and Figure 2 illustrate three-year time series of array yield (Yq) monthly average values, of 
the two PV technologies under study, from January 2018 to December 2020. A preliminary analysis of the 36 
data of each technology shows that the two arrays (pc-Si and a-Si) manifest a clear decreasing trend and a 
seasonal behavior. 

Figure 1 represents the pc-Si array yield variations over the 36 months from January 2018 to 
December 2020. The degradation of the pc-Si array yield is clearly shown by the decreasing trend of the curve. 
As for the variation due to the seasonality, the improvement of pc-Si performance is observed during the 
winter’s months. The high enhancement of pc-Si performance during the coldest months is due to the lowest 
module’s temperature, under the lower ambient temperature [21]. While during the warmest season the pc-Si 
efficiency is generally lower [22], [23]. The fluctuations of a-Si array yield (Yq) during the monitoring period 
are illustrated figure 2. As reported in many literature studies, the a-Si performance is influenced by Steabler 
Wronski effect [22] which explains the remarkable degradation of the array yield during the first year. 
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As reported by many literature studies, a-Si array reaches its peak of performance in high temperature condition 


due to the thermal annealing cycles; which helps recover some of performance loosed because of light induced 
degradation [24]. 


pce-Si Array Yield variations 


Array Yield YA 


0 5 10 15 20 25 30 35 40 
Months 


Figure 1. Array yield variations of pc-Si 


a-Si Array Yield variations 


Array Yield YA 


0 5 10 15 20 25 30 35 40 
Months 


Figure 2. Array yield variations of a-Si 


3.2. Degradation rates of array yield using statistical methods 
3.2.1. Array yields and classical seasonal decomposition method (CSD) 

The classical and seasonal decomposition is one of the most employed statistical techniques to 
estimate the performance loss. Using the additive model of this method, we can extract the trend, the seasonal 
and the irregular components from the initial time series of Ya the trend component helps get a fast idea of the 


performance degradation of the system [4]. Figures 3 and 4 illustrate the results found using the CSD technique 
under R: 
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Figure 3. Ya and CSD method for pc-Si 
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Figure 4. Ya and CSD method for a-Si 


The CSD method was applied to the time series of Ya, in order to quantify the degradation observed 
from the trend and described in section 3-1. As shown in Figures 3 and, according to this method the 
degradation rate for pc-Si system is found to be 0.44 %. While for a-Si the degradation rate is around 0.29%. 


3.2.2. Array yield and holt-winters model (HW) 

Although being less employed in PV degradation assessment, Holt-Winters method is also used to 
predict PV performance series [25].The HW model can be either additive or multiplicative, depending on the 
seasonal variations. As the seasonal fluctuations of PV systems are practically constant over the series, the 
additive model of Holt-winters algorithm is used in order to evaluate the degradation rates of the two PV 
systems [4]. The application of the linear regression to the HW series results in the determination of the slopes 
of the trends, which will be used according to (3) to calculate the performance loss rates. Figures 5 and 6 
represent results obtained after the application of Holt-Winter algorithm to the time series of the Array Yield 
(Ya) for the two technologies. 
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Figure 6. Ya and HW model for a-Si 


After the optimization of the Holt-Winters model under R, the parameters values of the model for pc- 
Si are o=0.11, B=0.49, y= 0.16. As for a-Si array, these parameters are o=0.03, B=1, y=1. The degradation rate 
DR (%) reaches for pc-Si system 0.99%. As for a-Si field, the YA depreciates at rate of 0.64%. 


3.2.3. Array yield and auto regressive integrated moving average method (ARIMA) 

ARIMA is also one of the most popular statistical methods, which is used; for PV systems; not only 
to assess the degradation of performance, but also to forecast the variation of performances in the future [26], 
[27]. Here, we are interested in the assessment of the degradation of the array yields of the two technologies 
under study. Again, the linear regression is called to evaluation the rates of the depreciation of performances. 


Figures 7 and 8 show results of application of the ARIMA method to the time series of Ya of the two fields 
under study. 
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Figure 7. Yq variations and ARIMA model for pc-Si 
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Figure 8. Y, variations and ARIMA model for a-Si 


After seeking the best ARIMA model that fits the times series under R, the parameters values for the model 
are ARIMA (2,1,0)(1,1,0)12 for pc-Si; while for a-Si field the parameters are ARIMA(0,1,2)(0,1,0)12. The 
degradation rate is around 0.48% for the pc-Si field; while for a-Si system it is about 0.32%. 


3.2.4. Array yield and seasonal-trend decomposition using LOESS (STL) 

The seasonal-Tend decomposition using LOESS is a powerful and flexible statistical method usually 
used in PV degradation assessment. The trend of the STL decomposition is explored in order to find out the 
slopes of the decreasing performance metrics of the two PV arrays. The application of the linear regression to 
the STL trends comes up with the estimation of the decline of array yields over time of the two PV systems. 
The Figures 9 and 10 exhibits the results of the application of STL method to the Yq time series of the two 
technologies. For pc-Si the degradation rate is found to be 0.39%, while for a-Si the Dr (%) is about 0,29 %. 
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Figure 9. Y, variations and trend of STL decomposition for pc-Si 
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Figure 10. Y, variations and trend of STL decomposition for a-Si 


3.3. Discussion of results 

Table 2 summarizes the annual degradation rates Dp (%) for the four methods, as well as the associated 
uncertainties and the coefficients of determination (R’). The different results of degradation rate Dr (%) 
acquired show that the degradation of PV performance depends on the statistical method used in addition to 
PV technology under consideration. Other parameters that affect the degradation rate are climate conditions 
[28], operating topologies, and the cumulative history of exposure to meteorological conditions [7]. 


Table 2. Degradation rates DR (%) of the two technologies under study 


Technology Parameter CSD Holt-Winters ARIMA STL 
Si DR(%) ~=—0.4440.05% 0.9940.19% 0.48+0.10%  0.39+40.04% 

aa R2 0.77 0.56 0.39 0.75 
a-Si DR (%) 0.29 40.03%  0.64+40.15%  0.3240.10%  0.29+40.02% 

R2 0.85 0.46 0.23 0.89 


In one side, the corresponding uncertainty to each degradation rate obtained from the application of 
the four methods give us the accuracy of each model. For both of pc-Si and a-Si arrays, the STL method 
provides the lowest Dr (%) with the best accuracy, 0.04% for pc-Si and 0.02% for a-Si. Here, STL model excels 
the CSD technique which comes up with low degradation rates with larger uncertainties: 0.05% for pc-Si and 
0.03% for a-Si. While the highest Dr (%) was given by HW model with the worst accuracy, 0.19% for pc-Si 
and 0.15% for a-Si. In the case of ARIMA method, it returns intermediate Dp (%) for the two PV arrays with 
intermediate uncertainty of 0.10% for both pc-Si and a-Si fields. 
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On the other side, values of coefficient of determination R* helps us to assess how well the linear 
regression model, used to extract the trend coefficient, fits the data generated by each technique. More the 
value of R* is close to ‘1’ more the degradation of the performance metric modeled is approximately linear. 
For the four techniques used in this work, the maximal values of R? are obtained by the application of STL 
method (0.75 for pc-Si and 0.89 for a-Si). Here again, STL outperforms CSD process which provides an R? of 
0.77 for pe-Si and 0.85 for a-Si. While intermediate and low values of the same coefficient R2 are presented 
by HW model (0.55 for pc-Si and 0.46 for a-Si) and ARIMA (0.39 for pce-Si and 0.23 for a-Si) models; 
respectively. This questions the linear approach conventionally adopted to deal with the degradation rate of 
different PV technologies [29]. 

As previously discussed, many factors influence the loss of performance of PV system. Table 3 (see 
Appendix) shows the degradation rates of different PV systems obtained by previous studies [10]-[12], [15], 
[30]+{33] using Several statistical methods, during various exposure periods and under dissimilar climate 
conditions. According to Képpen-Geiger climate classification [34], the weather in Laayoune City is 
considered to be a subtropical desert climate. A comparative study was carried out using results of nine previous 
works, which contain several PV plants. Only 17 systems are explored, which are composed of nine pce-Si and 
eight a-Si arrays. In one hand, various statistical methods are used to estimate the PV degradation rates. 
According to Table 3; LR/LLS and CSD are the most popular techniques to estimate the rate of performance 
losses Dr (%) (six studies for LR/LLS and four studies for CSD). On the other hand, high-performance losses 
can be found in composite climate (3.96% for a-Si), followed by desert and semi-arid climate (0.99% for pc- 
Si and 1.99% for a-Si). While low degradation rates are observed in semi-arid (0.2% for a-Si) and temperate 
climate (0.21% for pc-Si). If we do not consider the composite climate and look over the degradation rates of 
all parameter metrics; the loss of performance of a-Si fields passes by a minimum rate of 0.2% and a maximum 
of 1.99% with an average rate of 1.32%. As for pc-Si arrays the degradation rates have a minimum of 0.21% 
and a maximum of 0.99% with an average value of 0.61%. The comparative analysis conducted reveals that 
the degradation rates found in this work are in agreement with values mentioned in previous works. 


4. CONCLUSION 

In this paper, degradation assessment of the two silicon PV technologies (pc-Si and a-Si) over three 
years of outdoor-exposure and under climate conditions of Laayoune City has been conducted. The four 
statistical methods applied to Ya time series, using R software are: classical seasonal decomposition (CSD), 
holt-winters (HW), auto regressive integrated moving average (ARIMA) and Seasonal-trend decomposition 
using LOESS (STL). The main results and findings are recapitulated: 1) the two technologies present a seasonal 
behavior; ii) the annual degradation rate for pc-Si is found to be around 0.44%, 0.99 %, 0.48% and 0.39 %, 
using CSD, HW, ARIMA, and STL respectively with an average rate of 0.58%/Year; iii) while the array yield 
of a-Si degrades at rate of 0.29%, 0.64%, 0.32% and 0.29%, using CSD, HW, ARIMA, and STL respectively 
with an average rate of 0.39%/year; iv) Si array’s performance metric degrades less than that of pc-Si under 
the meteorological conditions of Laayoune; and v) the four statistical methods have shown that the two 
technologies exhibit an annual degradation rate less than 1%; with different uncertainties. The STL method 
excels the three other techniques in terms of accuracy (+0.04% for pc-Si and +0.02% for a-Si) while the highest 
uncertainty is obtained using the HW model (+0.19% for pc-Si and +0.15% for a-Si). vi) Results and findings 
of the present work are in agreement with results mentioned in previous works. 


APPENDIX 
Table 3. Degradation rates Dr(%) of PV systems from previous studies 

Location and Climate (Képpen - Technology of Pv Duration Parameter Methodology Dr Reference 
Classification) Module (Years) Metric (%) 
AIST Tsukuba pe-Si 3 PRr-25 LR -0.31 [30] 
Japan- pe-Si -0.51 
temperate climate pe-Si -0.55 

pe-Si -0.46 

pe-Si -0.21 

pe-Si -0.34 

sc-Si -0.76 

sc-Si -0.73 

sc-Si -0.96 

a-Si:H/c-Si -0.27 

a-Si -1.43 

a-Si -1.38 

CIGS 0.16 

CIGS 0.23 
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Table 3. Degradation rates Dr(%) of PV systems from previous studies (continue) 


Location and Climate (K6ppen - Technology of Pv Duration Parameter Methodology Dp Reference 
Classification) Module (Years) Metric (%) 
Tucson, Arizona, USA a-Si 3 Yf Bayesian statistical -0.20 [31] 
Semi-Arid analysis 
Saida, Algeria mc-Si 3 PR LLS -0.83 [32] 
Desert and semi arid (Multicristalin) P*M -0.74 
c-Si (mono- PR LLS -0.79 
cristallin) P*M -0.58 
HiT PR LLS -1.92 
P*M -1.53 
University of Cyprus in Nicosia,  pc-Si 3 PR CSD -0.99 [35 
Germany-Desert and semi arid ARIMA -0.92 
a-Si 3 PR CSD -1.87 
ARIMA -1.99 
Telangana, India mc-Si 4 PR LLS -0.27 [12 
Tropical CSD -0.32 
semi-arid HW -0.50 
STL -0.27 
Konya Anatolia, Turkey- pe-Si 2.5 P*M LLS -0.83 [10 
Semi arid HIT -0.67 
Jaén (Spain) a-Si 5.5 P*exp CSD “1.22 [il 
-Temperate YOY -1.32 
PRDC CSD -1.16 
YOY -1.40 
PRDC, STC CSD -1.25 
YOY -1.27 
India- a-si_l 3 PR LR -3.76 [30] 
Composite Climate Condition CSD -3.96 
LOESS -2.83 
a-si_2 LR -3.67 
CSD -3.36 
LOESS -2.11 
a-si_3 LR -3.73 
CSD -3.12 
LOESS -2.32 
HIT-1 LR -0.63 
CSD -0.48 
LOESS -1.33 
HIT-2 LR -0.26 
CSD -0.24 
LOESS -1.84 
HIT-3 LR -0.38 
CSD -0.34 
LOESS -1.87 
me-Si-1 LR -5.12 
CSD -5.04 
LOESS -2.70 
me-Si-2 LR -4.09 
CSD -3.84 
LOESS -2.40 
mc-Si-3 LR -4.29 
CSD -4.32 
LOESS -2.32 
Seville, (Spain)-Temporate mc-Si 22 Pmax LR -1.40 [15] 
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